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Abstract 

The natural neutron background influences the maximum achievable sensitivity in most deep underground nuclear, astroparticle and 
double-beta decay physics experiments. Reliable neutron flux numbers are an important ingredient in the design of the shielding of 
new large-scale experiments as well as in the analysis of experimental data. 

Using a portable setup of ^He counters we measured the thermal neutron flux at the Kimballton Underground Research Facility, 
the Soudan Underground Laboratory, on the 4100 ft and the 4850 ft levels of the Sanford Underground Research Facility, at the 
Waste Isolation Pilot Plant and at the Gran Sasso National Laboratory. Absolute neutron fluxes at these laboratories are presented. 

Keywords: Neutron flux. Underground laboratory, He-3 counters 


1. Introduction 

The major challenge of many modern physics experiments is 
the measurement of very small event rates (down to a few events 
per year, e.g. |[T|). Locating the equipment deep underground is 
a first step towards achieving this goal, since the rock overburden 
shields the experiment from cosmic rays. Generally the muon 
flux and the secondary gamma- and neutron fluxes are attenuated 
by a few orders of magnitude compared to the surface ||2l. 

A wide range of deep underground experiments are sensitive 
to background neutrons: In underground nuclear astrophysics 
stellar neutron sources need to be measured down to very low 
cross sections. Elastic scattering of neutrons can mimic sig¬ 
nals expected from WIMP interactions. Neutrinoless pp decay 
searches can be influenced by y rays emitted after neutron in¬ 
elastic scattering or capture and also by the decay of unstable 
nuclei produced through neutron capture. Therefore, the neu¬ 
tron background needs to be understood and sufficient shielding 
needs to be implemented to limit its impact on the experimental 
sensitivity. 

The underground neutron flux is mostly due to spontaneous 
fission of ^^*U in tfie cavity walls, {a, n) reactions induced by 
a-particles from the natural radioactivity of the underground 
environment and from the activity of the experimental setup 
itself OHni- The cosmic-ray induced neutron flux is two to 
three orders of magnitude lower than the radiogenic component 
Cl. Usually the laboratory neutron flux is simulated based on 
the composition of the rock and the concentration of radioactive 
isotopes iniiia. However, these simulations carry a large 
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degree of uncertainty and should be tested against measurements, 
if available. 

Geological composition of the rock and the varying contents 
of uranium and thorium in the underground environment as well 
as differences in the water content of the surrounding rocks 
cause variations in the background between underground labora¬ 
tories. In addition it has been found that even local differences 
in the composition of the rock can lead to background levels that 
vary by an order of magnitude between locations in the same 
laboratory (e.g.. Halls A and C in Gran Sasso) ||9l- 

Some data on measured neutron backgrounds are available, 
but a direct comparison is made difficult by the variety of detec¬ 
tion setups used and differences in the covered neutron energy 
range. In this work we present measurements of the thermal 
neutron fluxes at various underground locations using a portable 
setup of ^He detectors. Measurements were done at the Waste 
Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico; the 
Soudan Underground Laboratory in Minnesota; the Kimballton 
Underground Research Facility (KURF) in West Virginia; the 
Sanford Underground Research Facility (SURF) located in the 
Black Hills in South Dakota and the Laboratori Nazionali del 
Gran Sasso (LNGS) in Italy. 

2. Site background parameters 

The measurements presented here cover very different back¬ 
ground environments: WIPP is located in a salt formation, 
Soudan is a former iron mine, KURF is situated in an active 
limestone mine, SURF is located in a retired gold mine; the 
Gran Sasso laboratory is in a limestone formation. The prop¬ 
erties of each site relevant to the background conditions and 
measured values of the muon flux are listed in table Although 
WIPP is at a relatively shallow depth, the radiogenic radioac¬ 
tivity is very low due to the low-activity salt environment. U 
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Table 1: Properties of the visited underground sites. SURF values are stated for the 4850 foot level. See text for details on the WIPP U and Th concentrations. 



WIPP mi 

Soudan n2l|15| 

KijRF nm 

SURF iniini 

LNGS Him 

Environment 

Salt 

“Ely Greenstone” 

Limestone 

Poorman foundation 

Limestone 

Depth [m] 

655 

780 

500 

1500 

1400 

Equivalent depth [mwe] 

2000 

2090 

1400 

4300 

3800 

muon flux [lO^^ s^' cm^^] 

4.77 + 0.09 

2.0 + 0.2 

20 

0.044 + 0.001 

0.32 + 0.01 

[ppm] 

(0.48-1.49) -10-2 

0.17 


3.4 

6.8 (Hall A) 

232Th [ppm] 

(1.01 - 1.9) -10-2 

0.89 


7.1 

2.2 (Hall A) 


and Th concentrations displayed in the table were measured 
by inductively coupled plasma mass spectrometry on salt sam¬ 
ples. The low and high values from the range given in table 
were determined on a clear and a rocky sample showing a 
slight coloration, respectively. Additional measurements by 
gamma-spectroscopy are in agreement with these values and, 
by measuring the gamma-active daughter nuclides of the ura¬ 
nium and thorium decay chains, conhrm that they are in secular 
equilibrium. 

The thermal neutron flux at WIPP has been reported pre¬ 
viously as (1.3 + 0.3) ■ 10“^ cm“^ s“' lU. At LNGS various 
measurements have been done; none of them agree with each 
other, possibly due to a variation in the background depending 
on the location in the laboratory or due to unknown systematic 
uncertainties. The reported values are: (2.05 + 0.06) ■ 10“® cm“^ 
s-' 1181,(1.08 + 0.02)-10-® cm-2 s-> ||T^. and (5.4 + 1.3) ■ IQ-^ 
cm-^ s’ ESI. Only higher-energy neutron data are available at 
the other sites. 

3. Experimental setup 

3.1. Neutron detection with ^He proportional counters 



Figure 1: Typical spectrum from a ^He counter (1 channel « 1 keV). A neutron 
generates a signal betv^een channels ~200 and 800. Signals above and below 
this region are due to alpha particles and electronic noise, respectively. 

The ^He detectors used in our experiments consist of an 
A1 tube that is filled with ^He (and a small amount of CO 2 as 
quench gas) under a pressure of 10 bar. A wire at a -1-1400 V 
potential runs through the center of the A1 cylinder. ^He has a 


very high cross section for capmring thermal neutrons through 
the reaction ^He(n,p)^H (cr = 5330 barn, Q = 764 keV ED)- 
After a neutron has been capmred the proton (p) and the triton 
(T) deposit their kinetic energy (E^ = 573 keV, E7- = 191 keV) 
in the ^He gas. In case both reaction products are fully stopped 
in the sensitive volume a pulse with a height proportional to 764 
keV is generated. Due to the finite volume of the detectors there 
is a chance that one or both nuclei hit the A1 cylinder and only 
deposit a fraction of their energy in the detector (wall effect). 
This gives lead to the characteristic pulse height spectrum of 
a ^He proportional counter. Eig. [T] shows a typical spectrum 
taken with one of the ^He counters used in this work, with the 
two wall-effect peaks around channels 200 and 550. y-rays and 
electronic noise generate counts below channel 180 and can be 
clearly distinguished from the neutron signals. 

3.2. Intrinsic background 



Channel 


Figure 2: Low-gain spectrum. The thermal neutron peak is visible around 
channel 200. The counts at higher energies are due to the internal radioactivity 
of the counter itself. 

The counts visible at energies above the neutron peak are 
due to the internal radioactivity of the counters themselves, a- 
particles from the decay of the uranium and thorium present in 
the detector walls generate background signals covering the area 
of the thermal neutron peak (at 764 keV and below) and extend¬ 
ing up to 9 MeV Q. An additional background component is 
due to microdischarges near the central wire of the detectors 
II 22 II . The count rate in this range orresponds to the sensitivity 
limit for low-level neutron detection. The combined effect can 
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be seen in Fig. where the amplification of the detector signal 
has been lowered to cover a wider energy range. The thermal 
neutron full-energy peak lies at channel 200 with the internal 
detector background extending beyond it. This background is 
usually not of concern on the earth’s surface but at the very low 
neutron background conditions in an underground environment 
it becomes a major background component that can be stronger 
than the actual neutron signals. The average background rate 
integrated over the region of the neutron peak is about 10“^ s ' 
for the counters used here. Using the ratio of alphas in the neu¬ 
tron signal range to total alpha particles from the simulation and 
the area of the counter one obtains a total alpha activity of the 
counters of about 6 ■ 10“^cm“^s *. This value is in very good 
agreement with measurements of similar He-3 counters and an 
assay of commercially available aliuminum M- 

To model the a-induced background Geant4 simulations 
1231 were performed with a-particles being emitted from inside 
the aluminum container of the ^He counter. Contributions from 
decays in the anode material are negligible due to its much 
smaller surface area. Secular equilibrium in both the thorium 
and uranium chains was assumed and the particle energies were 
randomly chosen from the a-energies in the chains. Since the U 
and Th content of the containers is unknown various simulations 
with different Th/U ratios were performed; the resulting spectra 
are not significantly different from each other. A spectrum for 
Th/U = 1 is shown in Fig. the inset is a zoom into the energy 
region of the thermal neutron peak between 200 keV and 800 
keV. In the area of interest the energy deposition from alphas 
is featureless and the a-induced background can be reliably 
subtracted from the underground neutron data. 



Figure 3: Geant4 simulation of the pulse height spectrum due to cr-particles 
emitted from the walls of a ^ He counter. The region of the thermal neutron peak 
is marked by the dashed lines in the inset. 


3.3. Detector components and response function 

The experimental setup used for the underground ema sure - 
ments consisted of two ^He counter^as described in Sec. ITT] 


^The counters are 25.4 cm long (defined by the active part of the ^He inside 
the aluminum enclosure) and have a radius of 1.27 cm. 


GEANTTlsimulations resulted in a thermal neutron detection 
efficiency of 80 %. 

Details of the measurements at each site are described in the 
following section. In all cases the data were saved every hour so 
that the count rate could be checked for inconsistencies. 

4. Overview of the measurements the different sites 

4.1. KURF 

Data at KURF were taken for approximately one month. 
The detectors were deployed in the experimental hall on the 
14th level of the mine. At that depth the rock overburden is 
approximately 500 meters, equivalent to 1400 meters of water. 
Since the laboratory is located in an active limestone mine the 
study had to take into account nightly blast from the mining 
operations. No effect on the recorded count rate was observed. 

4.2. Soudan 

The experimental halls of the Soudan laboratory are located 
at a water equivalent depth of 2090 m (780 m below the surface) 
on the 27th level of the former Soudan Iron Mine. The detectors 
were set up on the second floor of the Soudan 2 / CDMS II 
cavern. Data were taken over a period of ~ 15 weeks. The 
detectors were calibrated with a ^^^Cf source. 

4.3. SURF 

Measurements at SURF were done on the 4100 foot and at 
two locations on the 4850 foot level. Three months worth of 
data were accumulated at 4100 feet. There the detector was 
placed inside an airlock along the main drift between the Yates 
and Ross shafts. The first measurement at 4850 feet was done 
for two months in an alcove just outside the entrance to the 
LUX and Majorana clean area. Then the detector was moved 
to the Majorana electroforming laboratory (“Temporary Clean 
Room”, TCR) and another 1.5 months of data were taken. After 
the underground measurements were completed the setup was 
brought to the surface and a short calibration run using the 
natural neutron background was performed. 

4.4. WIPP 

The neutron flux at WIPP was measured for approximately 
five months. The detectors were located in a connex in the “Q” 
alcove near the air intake shaft. WIPP lies at a depth of ~ 655 
m (ss 2000 m.w.e.), roughly in the center of a 600 m thick salt 
formation. A distance of at least 600 m separates the low-level 
transuranic waste storage areas from the experimental alcove, 
so that no adverse effect on the measured neutron flux is to be 
expected. 

^Version 9.6.2, with the high precision neutron data libraries G4NDL4.2, 
including thermal cross sections. 
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Figure 4: (color online) Spectra of the different measurements, (a) shows the results of the different SURF measurements in comparison with WIPR (b) compares the 
LNGS, Soudan and KURF data to WIPR The histograms in (b) were rebinned to easier differentiate them from each other. 


4.5. LNGS 

The setup was located in the upper connex of the LUNA II 
experiment, situated in the “interferometer tunnel” just outside 
of Hall A. The averaged rock overburden of LNGS is 1400 m 
(3800 m.w.e.) below the Gran Sasso mountain massif. The dom¬ 
inant rock is limestone. Before the measurement the detector 
calibration was checked using an ^"^'Am source. 

5. Analysis and results 

The thermal neutron flux (p can be calculated from the raw 
count rate Ctot. in the energy region of the neutron peak, the 
expected background count rate Chg, the detector efficiency tj 
and its surface area A = 213 cm^: 

Cfot Ghn 

cp ^ ' ( 1 ) 

rjA 

Figure shows the raw neutron data. To calculate the signal 
count due to neutrons the alpha counts from above the neutron 
peak were extrapolated to lower energies and subtracted from 
the raw counts since those resemble a mixture of neutron, alpha 
(and microdischarge) signals. The flux at WIPP is so low that 
even after 5 months of data taking no neutron signal is present. 
The spectrum therefore resembles pure background and can be 
used as a baseline level for comparison. 

The thermal neutron fluxes at each site are listed in Table 
1^ For WIPP and KURF it is only possible to determine upper 
limits (2 cr limits are given). Systematic uncertainties of +5% 
each for the efficiency simulation and the extrapolation of the 
alpha background into the signal region are included in the table. 

6. Discussion 

We presented thermal neutron flux measurements at the main 
underground laboratories in the United States and at LNGS, 


10 = 


5icr' 


Kb* ^ 


Soudan LNGS SURF 4100 SURF Davis SURF TCR KURF WIPP 

2090 mwe 3800 mwe 3600 mwe 4300 mwe 4300 mwe 1400 mwe 2000 mwe 

Site 


Figure 5: Neufron flux values or upper limits at the visited sites. The average 
flux on the earth’s surface is of the order of 10“^ s“* cm“^. 

presently the worlds largest undergound laboratory. Since identi¬ 
cal equipment was used to measure at each lab this data provides 
for the first time a systematic comparison of the relative neutron 
background levels at the different underground sites. Fluxes 
above a few 10“^ s * cm“^ have been determined with reason¬ 
able statistics with measurement times of a few weeks; the 
main limiting factor is the intrinsic radioactivity of the counters 
themselves, which can be clearly seen in the very long WIPP 
measurement. 

A comparison of the fluxes at each site is shown in Fig. 
The measured values lie in the expected ranges and show a re¬ 
duction in the flux of 3-4 orders of magnitude compared to the 
surface. We agree with the lowest reported value from LNGS 
II 20 II and reach a comparable uncertainty. The variation in flux 
between the three locations at SURF should not be seen as very 
surprising: there is no reason to expect that the neutron flux at 
a depth of approximately 1400 m is higher - or lower, for that 
matter - than on a mine level 250 m below. As stated in the 
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Table 2: Experimental results. The upper limits are 2 cr values, uncertainties are 1 cr. 



WIPP 

Soudan 

KURF 

LNGS 

Measurement time [10** s] 

12.7 

8.9 

2.3 

8.4 

Neutron flux [lO^^crn^^ s“'](+sfaf. + sys.) 

< 0.06 

0.7 + 0.08 + 0.05 

< 0.4 

0.32 + 0.09 + 0.02 


SURF (4100L) 

SURF (Davis) 

SURF (TCR) 


Measurement time [10*’ s] 

6.4 

4.9 

4.5 


Neutron flux [lO^^crn^^ s“'](+sfaf. + sys.) 

9.9+ 0.1 +0.7 

1.7+ 0.1 +0.1 

8.1 +0.1 +0.6 



introduction, due to the large reduction in cosmic rays at deep 
underground locations the neutron flux is entirely dominated by 
the radioactive components in the laboratory environments, e.g., 
the rocks and shotcrete near the experimental setup. Another fac¬ 
tor which can influence the local flux is the radon concentration 
in the vicinity, which is strongly dependent on the air circulation. 
In fact the highest neutron flux was measured at SURF’s 4100 
foot level, where the detectors were situated in a tightly enclosed 
volume with very little air flow. 

As is to be expected from a salt environment WIPP has by 
far the lowest neutron flux of the visited sites. Our upper limit is 
lower than the previously reported flux; aside from systematic 
differences this can be due to different concentrations of radioac¬ 
tive isotopes near the older measurement location. The local 
neutron flux is very sensitive to variations in the composition of 
the sutTounding environment. 
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